General response function for interacting quantum liquids 
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Linearizing the appropriate kinetic equation we derive general response functions including selfcon- 
sistent mean fields or density functionals and collisional dissipative contributions. The latter ones 
are considered in relaxation time approximation conserving successively different balance equations. 
The effect of collisions is represented by correlation functions which are possible to calculate with 
the help of the finite temperature Lindhard RPA expression. The presented results are applicable to 
finite temperature response of interacting quantum systems if the quasiparticle or mean field energy 
is parameterized within Skyrme - type of functionals including density, current and energy depen- 
dencies which can be considered alternatively as density functionals. By this way we allow to share 
correlations between density functional and collisional dissipative contributions appropriate for the 
special treatment. We present results for collective modes like the plasmon in plasma systems and 
the giant resonance in nuclei. The collisions lead in general to an enhanced damping of collective 
modes. If the collision frequency is close to the frequency of the collective mode, resonance occurs 
and the collective mode is enhanced showing a collisional narrowing. 



I. INTRODUCTION 

The response of matter to an external perturbation is 
the main source of knowledge about the matter itself. 
For instance, in plasma physics the polarization function 
II(<7, u>) is linked via the dielectric function e(q, u>) to the 
electrical conductivity a(q, ui) by 



e(q,cj) = l-V(q)U(q,ij) 



1 + ~a(q,u}). 



(1) 



In nuclear matter, e.g., the response function x = n/e al- 
lows one to study excitations and giant resonances which 
in turn yields information about the equation of state like 
the isothermal compressibility which is given by 
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or to calculate fluctuations and diffusion coefficients. 

Two lines of theoretical improvements of the response 
function can be found in print recently. The first one 
starts from TDHF equations and considers the response 
of nuclear matter described by a time reversal broken 
Skyrme interaction [Qj2|. The other line tries to improve 
the response by the inclusion of collisional correlations 
|}L|| and for multicomponent systems 0J^]. In this pa- 
per we want to combine both lines of improvements into 
one expression and derive therefore the response function 
from a kinetic equation including mean field (Skyrme) 
and collisional correlations. 

We consider here interacting matter which can be de- 
scribed by an energy functional ( mean field) £ originally 
introduced by Skyrme [p|JlO|] and the residual interaction. 
The latter one we condense in a collisional integral addi- 
tional to the TDHF equation. Then the response to an 



external perturbation will contain the effect of Skyrme 
mean field and additionally the effect of residual inter- 
action. While this schema and the results are of general 
interest for any interacting Fermi- or Bose system, we will 
only mention application examples from nuclear matter 
and plasma physics. For the latter one we might consider 
the energy functional as a parameterization of the self- 
energy in line with the philosophy of density functional 
theory. By this way we have the freedom to share the cor- 
relations between mean field like density functional pa- 
rameterizations and explicit collisional or dissipative like 
correlations which are condensed in a relaxation time. Of 
course, when deriving this parameterizations microscop- 
ically special care is required to avoid double counting of 
correlations. 

Specifically, we want to obtain the density, current and 
energy response Xj XJ, Xe of an interacting quantum sys- 
tem 
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(3) 

to the external perturbation y cxt provided the density, 
momentum and energy are conserved 

n(R, t) = (R)$j (R)) = s J2 /(P. *). 
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J(R,t) = ]T< 
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= «^p/(p,R, t), 

p 

E(R, t) = ]T(W(R, t)) = s J2 e(P. *)/(P, R> *)• (4) 

j p 

Here s is the spin-isospin degeneracy and we express the 
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observables in terms of the Wigner function /(p, R, t) 
which is related to the one-particle density operator p by 



/(p,R^5>^ R (p+f|p|p-f) 



(5) 



temperature T and a local mass motion momentum Q. 
These local quantities will be specified by the require- 
ment that the expectation values for density, momentum 
and energy are the same as the expectation values per- 
formed with /. 



and introduce the quasiparticle (Skyrme) energy e 



e(p,R, t) 
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(6) 



in the spirit of Landau theory s = 45. We will neglect 
the contributions of energy gain which arise from nonin- 
stantaneous collisions jnj . Here the energy functional or 
mean field (Skyrme) energy s is assumed to be parame- 
terized as [[12[ 
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A. Conservation laws 

From (||) we see that the conservation laws for density, 
momentum and energy are fulfilled if the corresponding 
expectation value of the collision side vanishes 



£p(/-/ Lc ') = ° 



(ii) 



Please note that the occurrence of current contributions 
~ J breaks explicitly the time invariance. These terms 
appear with the same coefficient t\ as the effective mass 
and energy contribution in order to ensure Galilean in- 
variance. The density dependence a^l deviating from 
the one arising by Skyrme three-body contact interaction 
has been introduced and compared with experiments in 



Taking this into account we can express the deviation of 
the observables <j) = l,p,e from equilibrium considering 
V = /-/o = /-/ le +/ le -/oas 
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II. DERIVATION OF GENERAL RESPONSE 
FUNCTION 

We start the derivation of the response from the quan- 
tum kinetic equation for the density operator in relax- 
ation time approximation ^ 



'p + i[£ + V cxt , P ] 



(9) 



where the relaxation is considered with respect to the lo- 
cal density operator /5 l c ' or the corresponding local equi- 
librium distribution function 

/.., PlI , i)=A (5t«p*i)j (10) 

with the (Fermi/Bose) distribution /q. This local equi- 
librium is given by a local chemical potential p, a local 



* The quasiclassical Landau equation follows from gradient 
expansion as 
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where we have performed Fourier transform t — > — ilo 
and r — > iq. In the last line we restrict to linear response 
of @. We assume for simplicity a homogeneous equi- 
librium /o(eo(p)) such that only the deviations 5<p(r,t) 
and Sf(p,r,t) are space dependent. This is no principle 
restriction but otherwise a lot of later algebraic expres- 
sions would take the form of integral equations. Further 
<5Q(q, t) = q<5Q(q, t) is employed. With the abbreviation 
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and the correlation function 

/o(p + f)-/o(p-f) 
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(14) 



we can write the deviation of the observables from equi- 
librium according to (Il2|) explicitly 



Sn = — Suai — SQbi — STc\, 
SJ q = q<5J = -<5^a qp - <5Q& qp 
<5_E = —Sfia e — SQb e — STc e . 
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Instead of the vector equation for the current J we 
consider the projection onto the direction of q. This 
simplifies matters as long as we have no active media 
and J| |Q. 
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The matrices are 
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and the abbreviations are introduced 
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/ giVa+g^Vi 3 pq Vi giV 2 \ 

I ffpqVb + g p 2 p qV 4 ^pq) 2 ^ ffpq^J 
\ geVo + +Ve^4 ffepq^l ffc^2 / 



(20) 



B. Response from kinetic equation 

To derive the response function we will obtain a second 
equation set from linearizing the kinetic equation (|^) and 
the corresponding balance equations. Fourier transform 
t — > — ioj and r — > iq the equation (H) can be linearized 
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In ( pif ) the mean field contributions Vi will lead just to 
selfconsistency. The coefficients V, it-selves are linked to 
the parameterization urn as 
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in terms of the correlation functio n (|l4| ). The required 
solution is obtained from (p~8h and (|19[) as 



6j q I = (x- v-ka-^ -1 

<5£ 




(22) 



from which one can read off the response functions (|^). 
This is the main result of the paper which represents 
the density, momentum and energy response including 
nonlinear mean fields and collisions with the fulfillment 
of density -, momentum -, and energy - conservation. 



C. Alternative expressions 

Before we continue to consider special cases we like to 
express the solution (E2J) in a slightly more familiar form. 



We can solve (jig ) for Sf and perform momentum in- 
tegrations to obtain the observables Sn, 5J q , SE. This 
leads to the following closed equation system 
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SJ q 
SE 
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1. Response in terms of mean- field response 

First we assume that we have solved the response with- 
out collisions B = which would obey the equation 



(1 - V)£o = Qv e 



(23) 



together with the set ([L5h 



where £ = {Sn, SJ q , SE}, v ext = {V oxt , 0, 0} and 
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Than the response matrix (|3|) without collisions but self- 
consistent mean field reads 



X MF (co) = (1 - v)- 1 g{u). 



(25) 



The missing part of the full solution of ( 18|19| ) including 
the collisions are given by £ = £o + C where we have for £ 



(1 - V - BA- V K = BA- 1 ^- 
Some algebra leads to the final response 
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2. Response in terms of polarization function 

The opposite case is the usual way where we first solve 
the equation without selfconsistency by the mean field. 
This leads to the polarization function V = {II, IIj, II^;} 
which we use to represent the response function which in- 
cludes selfconsistency. Without mean field we have from 



(1-BA- 1 )^ = Qv exi 
leading to the polarization function 

V{uo) = (I — BA~ 1 )~ l G{tjj - 



(28) 
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The response function can be represented analogously to 



X(u) = V{tu) { 1 - Q- X {u + -) VV(ui) 



(30) 



The generalization of the usual form \ = 11/ (I — VII) for 
simple mean fields can be recognized. 



III. CALCULATION OF RESPONSE FUNCTIONS 

In the following we consider some frequently occurring 
situations. Therefore we assume only quadratic disper- 
sions e = p 2 /2m in the correlation function. To consider 
the full quasiparticle (Skyrme) energy would correspond 
to the selfconsistent quasiparticle RPA which we do not 
want to consider here. The schema how to include this is 



clear after the following considerations. We understand 
in the following the mass as effective mass given by t\ in 

The different occurring correlation functions (tL4j) can 
be written in terms of moments of the usual Lindhard 
polarization function Hq 
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For practical and numerical calculations we can rewrite 
the II„ by polynomial division into 
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where the IT are the projected moments perpendicular 
to q and read 



IT 



pq x 2 /o(p + f) -/o(p- 



(2^ (P "^ q) = ~ 



2m J dn'n 



n 4 = 



2mTn 
dp 



, PQ , 4 /o(P + f)-/o(p-f) 

(p - -o-q) i 



(2ir)° q 



uj — iO 



2(2m) 2 J dp! J dn"U 



8m 2 T 2 H . 



(34) 



The corresponding last identities are valid only for nondc- 
generate, Maxwellian, distributions with temperature T. 
The general form of polarization functions is presented as 
an integral over the chemical potential [i of the Lindhard 
polarization n . This is applicable also to the degen- 
erate case. In the following we will discuss successively 
further involved results; first for nondegenerate plasmas 
and then for degenerate nuclear matter. 
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A. Polarization with collisions: inclusion of density 
and energy conservation 



B. Polarization with collisions: inclusion of density 
and momentum conservation 



Now we concentrate on the response function without 
mean field and consider only the collisions within density 
and energy conservation. Then the matrices A and B 
reduces to 2x2 matrices. The calculation of (E8h leads to 
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(36) 



With the help of (|32j-|34|) this can be further worked out 
in terms of II n but leads not to a more transparent form. 



Let us note that the first term in (35) represents just the 
result if we would have considered only density conserva- 
tion known as the Mermin polarization function Q 



n n (w) = (1 - iter) 



ffi(^+7)ffi(°) 



(37) 



Of course the limit of vanishing collisions r — > oo en- 
sures that the Lindhard result Ho(lo) appears since 
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FIG. 1. The dielectric function e = 1 — VH for a one compo- 
nent plasma system in different approximations. The inverse 
Debye length is a 2 = 4 ™ e and the energy co is scaled in 
plasma frequencies ujp — K 2 T/rn. The relaxation time is cho- 
sen arbitrarily as r = 5/u> p and the wave vector q = 0.384k. 
The upper most panel shows the excitation function. 



Next we consider the special case that the density and 
momentum are conserved. Then the matrices A and B 
reduces again to 2x2 matrices and the calculation of ( |2"§| ) 
leads to 



n n ' J ( w ) = (i — iut) 



, , J$ 9i(0). (39) 

h _ (W*+-)) 



We have used the fact that according to (|32|) g pq (0) = 
.9epq_(0) = and Ha is defined as in (36). With the help 



of (|lj) we obtain by this way a slightly modified Mermin 
dielectric function (37). 



C. Polarization with collisions: inclusion of density, 
momentum and energy conservation 

Considering all three conservation laws the result from 
©is 
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This result together with the former special cases 
([il-ill'H) are compared in figure]!]. One sees that the 
first approximation of Mermin ( |37|) is almost identical 
with the result (|35| ) where density and energy is con- 
served. The inclusion of density and momentum conser- 
vation (|3^) brings the curves towards the Lindhard result 
without collisions compared with the inclusion of den- 
sity conservation only. Finally, the complete result with 
the inclusion of density, momentum and energy conser- 
vation (^) changes the results again in the direction of 
the result for density and energy conservation but less 
pronounced. This qualitative behavior of the different 
approximations are observed for other relaxation times 
too. 
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The effect of relaxation times within the complete re- 
sult (|^) is seen in hgure ||. One recognizes that with de- 
creasing relaxation time or increasing collision frequency 
the plasmon peak is shifted towards smaller energies. For 
collision frequencies around the inverse plasma frequency 
there occurs a resonance seen in the real part of the di- 
electric function (middle part of |^). This translates into 
an enhanced single particle damping (lower panel) and 
the system becomes optical thick. At the same time 
the collective mode, the plasma frequency, becomes en- 
hanced. One can consider this as an effect of transferring 
collisional energy into collective motion. We have here 
a coherent superposition between collision frequency and 
collective frequency resulting into an enhancement of col- 
lective motion. Compared with the general effect of col- 
lisions to increase the damping of collective motion, see 
figure [I], this is the inverse effect which narrows the col- 
lective peak again. 
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FIG. 2. The dielectric function e = 1 — VTI for a one compo- 
nent plasma system with the inclusion of density, momentum 
and energy conservation (^) for different relaxation times. 
The parameter are the same as in figure |l[ 



D. Response with collisions: simple mean field 

For the case of simple but density - dependent mean 
fields Vo 7^ and V\ = V% = V4 = we obtain for the 
density response from (pCf) 



n 



x 



1-VnII 



(43) 



Here II is the polarization function without mean field 
but with collisions. Dependent on the choice one may 
use ©, © or @ for the latter one. 



Since the imaginary part of the response function M3) 
is related to the photoabsorbtion yield on nuclei we like 
to apply the different approximations ([35]), ( p9| ) or fllc| ) 
for nuclear isovector oscillations. We use first a simpli- 
fied Skyrme parameterization [Q for Vb according to 
(Mq) and a relaxation time r within a Fermi liquid model 
|l5f,^6|. From Fig. |^ we see the same qualitative behav- 
ior of the different approximations as found in Fig. 
The difference between the density (density-momentum) 
and density-energy (density-energy- momentum) result is 
very small (inlays of Fig. |^) . An increase of temperature 
leads to larger damping of all approximations (bottom of 
Fig. H) . While the density or density-energy result leads 
to a pronounced damping of the giant resonance the in- 
clusion of momentum conservation diminishes this effect 
again towards the free result. 

Let us note that the energy-weighted sum rule (EWSR) 



du u> Imx 



2m 



"0 



(44) 



is fulfilled numerically for all approximations, however, 
the convergence is very bad for response including den- 
sity or density and energy conservation. The inclusion 
of momentum conservation in turn improves the conver- 
gence of the sum rule appreciable. 
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Mermin (n) 
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10 20 

(0 [MeV] 

FIG. 3. The imaginary part of the response function Imx 
for nuclear giant dipole resonances in different approximations 
(comp. Fig. |l|). The wave vector q = 0.23 fin -1 (« 0.17k F ) 
corresponds to the inverse diameter of the nucleus 208 Pb ac- 
cording to The inlays show an enlarged view of the 
difference between the density (n), density- energy (n-E) and 
density-momentum (n-p), density- momentum-energy (n-p-E) 
approximation, respectively. 



E. Response without collisions: Skyrme mean field 

First we consider the case that we have only Skyrme 
- mean fields. Then the matrix equation (E3h is solved 
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with the result for the density response function 
Xmf(w, Vi) = 
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For isovectorial oscillations one has 
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where the Vf are representing the Skyrme parameteriza- 
tion in nuclear matter Q 

V ' = -to(xo + \)-j(x 3 + ±)r% 
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If we use the definitions of which are related to ours 
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we obtain the result of [jlj which was slightly misprinted 
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for nuclear matter density no and s = 4. 

F. Response with collisions: inclusion of density 
conservation and Skyrme mean field 

Now we derive the combined result from the mean field 



response (45) and collisions. We restrict first only on 



density balance conservation [first term of (|3g)] to get 
from (flS 
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and in 
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Therefore we can solve (|30|) and obtain similar to (|45j 
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(55) 



We should note that the mean field potential Vi arising 
from current and effective mass is of a different level of 
approximation than the collisional contribution which we 
restricted in ( |53] ) to the density conservation. This incon- 
sistency is visible in the final result for Vi in (|5^) where 
the limit of vanishing collisions r — > 00 does not agree 
with the limit of finite collisions since g pq (0) = 0, see 
(|32|). Therefore Vi = V is proposed to ensure the limit 
of vanishing collisions. 
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FIG. 4. The imaginary part of the response function Imx 
for nuclear giant dipole resonances (see Fig. in different 
approximations. Considering the Skyrme interactions SGII 
MfiM we compare the full mean field response with collisional 
correlations in density (dotted line) and density, momentum 
and energy approximation (dashed line), respectively, with 
the collision-free mean field response (solid line) for two dif- 
ferent temperatures. 



G. Response with collisions: inclusion of density, 
momentum and energy conservation and Skyrme 
mean field 

The Skyrme response can be given also for the other 
cases including energy and momentum conservation. 
However this leads not to a more transparent form than 



the general matrix structure (30). Considering the stan- 
dard effective Skyrme interaction SGII H| in figure ^ 
we compare the complete result ( p2[) including energy, 
momentum and density conservation (dashed line) with 
the result including only density conservation (|54| ) (dot- 
ted line). In the latter one we have proposed for V\ the 
collision-free value V\ in order to ensure the correct lim- 
iting case. 

We find again the same behavior for the different ap- 
proximations as in the case of the simplified mean field 
(Fig. |^) . The inclusion of collisions leads to an enhanced 
damping and a shift of the collective peak towards smaller 
energies. This effect of collisions is less pronounced by the 
full mean field result including density, momentum and 
energy conservation (solid line). The increase of temper- 
ature leads to a broadening of the resonance structure 
(lower part of Fig. |^) in any case. 

Again we have checked our results to satisfy the EWSR 



du> U) Imx = 7, — n o(l + «) 
2m 



with the enhancement factor re 



[t 1 (2 + x 1 )+t 2 {2 + x 2 )]n . 



(56) 



(57) 



IV. SUMMARY 

In this paper we have derived the unifying response 
function including nonlinear mean fields (Skyrme) and 
collisional correlations. Within this approach one can 
share correlations between an energy functional of mean 
field (Skyrme) parameterization and explicit dissipative 
correlations condensed in the relaxation time. This al- 
lows to treat also dissipative effects in density functional 
approaches. We see that the known limiting cases are re- 
produced neglecting either collisions or mean fields. Spe- 
cial transparent cases of the unifying response are dis- 
cussed. 

For a nondegenerate plasma numerical results are pre- 
sented. The first order correction given by the Mermin 
response, incorporating only the density balance, is sim- 
ilar to the approximation where density and energy are 
conserved. The plasmon peak is shifted towards smaller 
frequencies. This is accompanied by an enhanced damp- 
ing. The incorporation of momentum balance diminishes 
this effect of collisions. 

We observe that an enhancement of the collective 
mode occurs for collision frequencies near the collective 
(plasma) frequency. This is the inverse effect of damping 
due to collisions in that the collisions become resonant 
and the collective mode is enhanced. We consider this as 
collisional narrowing. Since the momentum conservation 
is responsible for that effect we suggest that the physical 
origin is the same as sometimes discussed with motional 
narrowing. We like to stress that this narrowing is ob- 
served relative to the broadened mode due to collisions 
and did not reach the collision - free value. Consequently 
we have collisional damping every time but near the res- 
onant situation this collisional damping is diminished. 

Similar behavior is found for the case of nuclear mat- 
ter, where the collective mode is the giant resonance. For 
isovectorial giant resonances we checked the extended en- 
ergy weighted sum rules and find excellent completion. 
The response due to nonlocal mean fields is derived in- 
cluding the effect of collisional correlations. 



Here re occurs as a consequence of the momentum- 
dependent terms in the Skyrme interaction and is defined 
as the deviation from the Thomas-Reiche-Kuhn sum rule 
in the case of isovector giant dipole resonance [ ^9pc| ] . 

The result ( |45| ) which contains the full Skyrme mean 
field but no collisions is in excellent agreement with (56). 
The approximation ( p4^ ) including collisions but only den- 
sity conservation conserves the sum rule only w 75% 
which is a consequence of the inconsistency of (|55|) . The 
inclusion of density, momentum and energy conservation 
(p2|) conserves the sum rule (pw again completely. 
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